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Abstract 



The gluino-induced contributions to the decay b ^ sj are investigated in 
supersymmetric frameworks with generic sources of flavour violation. It is 
shown that, when QCD corrections are taken into account, the relevant op- 
erator basis of the Standard Model effective Hamiltonian gets enlarged to 
contain: i) magnetic and chromomagnetic operators with a factor of Os and 
weighted by a quark mass rub or mc, ii) magnetic and chromomagnetic op- 
erators of lower dimensionality, also containing a^; iii) four-quark operators 
weighted by a factor a^. Numerical results are given, showing the effects of 
the leading order QCD corrections on the inclusive branching ratio for b — > sj. 
Constraints on supersymmetric sources of flavour violation are derived. 



October 1999 



1 



I. INTRODUCTION 



Processes involving Flavour Changing Neutral Currents (FCNCs) provide invaluable 
guidelines for supersymmetric model building. The experimental measurements of the rates 
for these processes, or the upper limits set on them, impose in general a reduction of the 
large number and size of parameters in the soft supersymmetry-breaking terms present in 
these models. Among these processes, those involving transitions between first- and second- 
generation quarks, namely FCNC processes in the K system, are considered as the most 
formidable tools to shape viable supersymmetric flavour models. Moreover, the tight ex- 
perimental bounds on some flavour-diagonal transitions, such as the electric dipole moment 
of the electron and of the neutron, as well as — 2, help constraining soft terms inducing 
chirality violations. 

Several supersymmetric models have so far emerged, with specific solutions to the chiral- 
flavour problem. Among them are two classes of models in which the dynamics of flavour sets 
in above the supersymmetry breaking scale and in which the subsequent flavour problem is 
killed by the mechanisms of communicating supersymmetry breaking to the experimentally 
accessible sector. They are known as mSUGRA, i.e. minimal supersymmetric standard 
models in which supergravity is the mediator between the supersymmetry-breaking sector 
and the visible sector and gauge-mediated supersymmetry-breaking models (GMSBs) 0, 
in which the communication between the two sectors is realized by gauge interactions. In 
other classes of models, particular flavour symmetries are introduced, which link quarks 
and squarks: models in which an ahgnment of squarks and quarks is assumed 0, and 
models in which the solution to the flavour problem is obtained by advocating heavy first- 
and second-generation squarks In the latter, the splitting between squarks of first 

and second generation and those belonging to the third generation relies on a f/(2) flavour 
symmetry 

Neutral flavour transitions involving third-generation quarks do not yet pose serious 
threats to these models. One exception comes from the decay h — >■ 57, the least rare 
flavour- and chirality-violating process in the B system. It has been detected, but the 
precision of the experimental measurement of its rate is not very high at the moment. 
Nevertheless, this measurement already has the effect of carving out some regions in the 
space of free parameters of most of the models in the above classes (see for example [§; 
for a recent analysis, see P] and references therein). They also drastically constrain several 
somewhat tuned realizations of models in these classes |TO|JTT[] . Once the precision in the 
experimental measurement has increased, this decay will undoubtedly gain efficiency in 
selecting the viable regions of the parameter space in the above classes of models and it 
may help discriminating among the models by then proposed. It is therefore important 
to get ready reliable calculations of this decay rate, i.e. calculations in which theoretical 
uncertainties are reduced as much as possible, and which are general enough to be applied 
to generic supersymmetric models. 

The experimental situation is, at present, as follows. The ALEPH Collaboration at LEP 
reports a value of the inclusive decay B of [|T2[: 



BR(5 ^ X,7) = (3.11 ± 0.80 ± 0.72) x 10"^ (1) 
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from a sample of b hadrons at the Z resonance. The CLEO Collaboration at CESR has a 



statistically and systematically more precise result, based on 3.3 x 10 BB events [|13 
BR{B X,7) = (3.15 ± 0.35 ± 0.32 ± 0.26) x 10~^ , 



(2) 



but quotes a still very large interval 

2 X 10"^ < BR{B X,7) < 4.5 x 10'^ 



(3) 



as the range of acceptable values of branching ratios. 



characterized by its large QCD contributions, 
1^ , is known with high accuracy in 



Theoretically, the rate for this decay 
practically as large as the purely electroweak ones 
the Standard Model (SM). It has been calculated up to the next-to-leading order (NLO) 
in QCD, using the formalism of effective Hamiltonians [|15|. Results for LO and NLO 
calculations and for power corrections can be found in [|T^-|T8|, [|T^-|23|, and [Q, respectively. 



The resulting theoretical accuracy is rather astonishing: the inclusion of the NLO QCD 
corrections reduces the large scale dependences that are present at LO (±25%) to a mere 
per cent uncertainty, once the value of the parameters to be input in this calculation is fixed. 
This accuracy, however, is obtained through large and accidental numerical cancellations 
among different contributions to the NLO corrections and a subsequent cancellation of scale 
dependences p3| , |25|l . The same accuracy, indeed, is not obtained for the NLO calculation of 
the rate BR(i? —>■ X^j) in simple extensions of the SM, such as models that differ from the 
SM by the addition of two or more doublets to the Higgs sector |2^ . 

The calculation of BR(i? Xg'y) within supersymmetric models is still far from this 
level of sophistication. There are several contributions to the amplitude of this decay, usually 
identified by the particles exchanged in the loop. Besides the Vr~-t-quark and i/~-t-quark 
contributions, there are also the chargino, gluino and neutralino contributions, respectively 
mediated by the exchange of chargino-up-squarks, gluino-down-squarks and neutralino- 
down-squarks. All these contributions were calculated in Ref. within mSUGRA; their 
analytic expressions apply naturally to GMSB models also. The inclusion of QCD corrections 
needed for the calculation of the rate, was assumed in p6[| to follow the SM pattern. No 
dedicated study of this decay exists for the supersymmetric models mentioned above with 
specific flavour symmetries. A calculation of BR(i? Xg^) induced solely by the gluino 



contribution has been performed in p7| , |28[| for a generic supersymmetric model, but no QCD 
corrections were included. 



A NLO analysis of BR(i? — > Xs'y) was recently performed for a specific supersym- 
metric case (the corresponding NLO matching conditions are also given in 13^]). This is 
valid in a class of models where the only source of flavour violation at the electroweak scale 
is that of the SM, encoded in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. It applies 
to mSUGRA and GMSB models (in which the same features are assumed/obtained at the 
messenger scale) only when the amount of flavour violation, generated radiatively between 
the supersymmetry-breaking scale and the electroweak scale, can be neglected with respect 
to that induced by the CKM matrix. It applies, therefore, to the case in which only the 
lightest stop eigenstate contributes to the chargino-mediated loop and all other squarks and 
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gluino are heavy enough to be decoupled at the electroweak scale. It cannot be used in 
particular directions of parameter space of the above listed models in which quantum effects 



induce a gluino contribution as large as the chargino or the SM contribution |TT| , |32| . 
Nor can it be used as a model-discriminator tool, able to constrain the potentially large 
sources of flavour violation typical of generic supersymmetric models. 

Among these, flavour-violating scalar mass terms and trilinear terms induce a flavour 
non-diagonal vertex gluino-quark-squark. This is generically assumed to provide the dom- 
inant contributions to quark- flavour transitions thanks to its large coupling Qs- Therefore, 
it is often taken as the only contribution to these transitions |^^, and in particular to the 
6 — >• S7 decay, when attempting to obtain order-of-magnitude upper bounds on flavour- 
violating terms in the scalar potential p7| , |25| . Once the constraints coming from experimen- 



tal measurements are imposed, however, the gluino contribution is reduced to values such 
that the SM and the other supersymmetric contributions can no longer be neglected. Any 
LO and NLO calculation of the — 57 rate in generic supersymmetric models should then 
include all possible contributions. 

The gluino contribution presents some peculiar features, related to the implementation 
of QCD corrections, that have not been detected so far. As already mentioned, the decay 
6 — s> S7 involves a quark-flavour violation as well as chirality violation. The first is directly 
related to the flavour violation in the virtual sfermions exchanged in the loop. The second can 
be obtained as in the SM, through a chirality flip in the external 6-quark, and it is signalled 
by its mass m^. It can also be induced by sfermion mass terms originating from trilinear 
soft supersymmetry-breaking terms. These mass terms differ from fermionic mass terms by 
two units of i?-charge under a U{1)r symmetry. The correct i?-charge for this b-s transition 
is then restored through the insertion of the gluino mass m~ in the gluino propagator. 
The two different mechanisms producing chirality violation are well known. They give 
rise to operators of different dimensionality when generating the effective Hamiltonian used 
to include QCD corrections to the 6 — 57 decay. Indeed, m~, the mass of one of the 
heavy fields exchanged in the loop, is naturally incorporated in the Wilson coefficient of 
the corresponding magnetic operator, which is now of dimension five {egl {sa'^'^Pub) F^iy). 
On the contrary, Wb, the running mass of one light field, with a full dynamics below the 
matching scale, is naturally included in the definition of a magnetic operator, which is of 
dimension six {e g1mi,{sa^'' PrI)) F^y). 

Moreover, the presence of the strong coupling as in the gluino contribution immediately 
sparks off the question of whether this coupling should be included in the definition of the 
gluino-induced operators or in the corresponding Wilson coefficients. Both choices are, in 
principle, acceptable. It can be observed, however, as will be discussed in Sec. |I|, that 
the first option does not require a modification of the program of implementation of QCD 
corrections established in the SM case. In particular, the anomalous dimension matrix 
starts at order as and is used up to order as (a^) in a LO (NLO) calculation. The inclusion 
of the as coupling in the operators imposes a necessary distinction of the dimension six 
gluino-induced magnetic operators e g1miy{sa^^ Pr})) F^j^y from the SM magnetic operator 
e/167r^mfe [sa^'^Pub) F^y. As it will be seen in Sec. a set of new four-fermion operators, 
induced by gluino exchanges, is also needed. 
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These features single out the gluino contribution to the decay 6 — 57 as one that neces- 
sarily requires a dedicated study of the implementation of QCD corrections already at the LO 
in QCD, before including chargino and neutralino contributions and higher-order QCD cor- 
rections. In Sec. H, the list of operators induced by gluino-mediated loops is given together 
with the list of those needed for the SM contribution. The number of operators depends 
on the sources of flavour violation that are present in the particular supersymmetric model 
considered. In the attempt to reach the level of generality advocated above, no restriction is 
made on the possible sources of flavour violation in the sfermion sector. These are surveyed 
in Sec. |T|. Also shown is the direct connection between flavour-violating sources and opera- 
tors generated, emphasizing the differences between the analysis in a generic supersymmetric 
model and the typical mSUGRA-inspired analyses. The Wilson coefficients at the matching 
scale for the Hamiltonian generated by gluino contributions are given in Sec. They are 
calculated using the mass-eigenstate formalism, the most appropriate to deal with different 
off-diagonal terms in the sfermion mass matrix squared, of a priori unknown size. These 
coefficients evolve down to the low-scale fib independently of the usual SM coefficients, since 
there is no mixing between SM and gluino-induced operators. The anomalous-dimension 
matrix governing this evolution at the LO in QCD and the resulting analytic expressions 
for the low-scale Wilson coefficients is given in Sec. |V|. In Sec. |VJ an expression for the LO 
rate BR{B — >• Xg'j), due to the SM and the gluino-induced Wilson coefficients, is derived. 
Numerical evaluations of the branching ratio are shown in Sec. |VII| , when only one or at 
most two off-diagonal elements in the down-squark mass matrix squared are non- vanishing. 
As already mentioned, the decay 6 — >^ 57 can be realistically used as a tool to select viable 
supersymmetric flavour models only when all contributions to BR.{B Xg'y) are included. 
The numerical evaluations of Sec. |VII| , therefore, have only the purpose of illustrating the 
effect of the LO QCD corrections, as well as the interplay between SM and gluino contri- 
butions to the branching ratio. Strictly speaking, they give results that are valid only in 
particular directions of the parameter space of generic supersymmetric models, and provide, 
in general, some intermediate results of an ongoing, more complete analysis. 

II. ORDERING THE QCD PERTURBATIVE EXPANSION AND THE 

EFFECTIVE HAMILTONIAN 

In the SM, rare 5-meson decays are induced by loops in which W bosons and up-type 
quarks propagate. The most important corrections are due to exchanges of light particles, 
gluons and light quarks, which give rise to powers of the large logarithmic factor L = 
log(mg/m^). 

The decay amplitude for b s'j obtains large logarithms L only from loops with gluons. 
This implies at least one factor of for each large logarithm. Since the two scales rrib and 
Mw are far apart, L is a large number and these terms need to be resummed: powers of 
agL are resummed at the LO, terms of the form ag (agL)^ are obtained at the NLO. Thus, 
the corrections to the decay amplitude are classifled according to: 

(LO): G^(a.L)^, (AT = 0,1,...) 
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(NLO): GpasiasL) 



N 



where is the Fermi constant. 

The resummation of these corrections is usually achieved by making use of the formalism 
of effective Hamiltonians, combined with renormalization group techniques. The needed 
effective Hamiltonian is obtained by integrating out the heavy degrees of freedom, i.e. the 
top-quark and the W boson. It is usually expressed as 

AC 

nZf = -^v,,vz y: , (4) 

where Vtb and Vts are elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The 
Wilson coefficients Cj contain all dependence on the heavy degrees of freedom, whereas the 
operators Oi depend on light fields only. The operators relevant to radiative B decays can 
be divided into two classes: 



• current-current operators and gluonic penguin operators |T^: 

04 = is^^.T-PLb)j:giqrTaq), 

= {si,J.%PLb)j:g{qrrYq), 
Oe = {-si^.lulpT-PLh)Y.,{qYYYTaq) , 

where T" (a = 1,8) are SU{'i) colour generators; 

• magnetic operators, with chirality violation signalled by the presence of the 6-quark mass: 

^7 = ^^m,{^^){m^''PRh)F^,, 

(6) 

Os =^,m,ifx)isa^-'T^Pnb)G';,,, 

where Qs and e are the strong and electromagnetic coupling constants. Both sets of operators, 
those in and in (|^) are of dimension six. 

It is by now well known that a consistent calculation for b 57 at LO (or NLO) precision 
requires three steps: 

1) a matching calculation of the full standard model theory with the effective theory at 
the scale n = fj,w to order q;° (or al) for Wilson coefficients, where fiw denotes a scale 
of order Mw or rrit; 

2) a. renormalization group treatment of the Wilson coefficients using the anomalous- 
dimension matrix to order a] (or a^); 
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3) a calculation of the operator matrix elements at the scale fi = fib to order a° (or a]), 
where fib denotes a scale of order mb- 

That matters can be somewhat different is illustrated by the decay b — sii. The 
effective Hamiltonian @ contains in this case two additional operators: 



e2 



^9 =T7—Asi,PLb){W 



^2 



(7) 



010 = TTT^ {si,PLh) {h^ld) 



It turns out that in this case, the operator O2 mixes into O9 at one loop: the pair cc in O2 
can be closed to form a loop, and an off-shell photon producing a pair can be radiated 
from a quark line. The first large logarithm L = log(m^/M^) arises without the exchange 
of gluons. This possibility has no correspondence in the b ^ s'y case. Consequently, the 
decay amplitude is ordered according to GpL (asL)^ at the LO in QCD and GpctsL^asL)^ 
at the NLO. To achieve technically the resummation of these terms, it is convenient to 
redefine magnetic, chromomagnetic and lepton-pair operators Ot, Og, Og, and Oiq and the 
corresponding coefficients as follows 



0-- = i^O,, cr = S^a (^ = 7,...,10). (8) 



This redefinition allows us to proceed according to the above three steps when calculating 
the amplitude of the decay b ^ s££ |^ . In particular, the one- loop mixing of the operator 



O2 with the operator Og'^'^ appears formally at 0{as)- 

In supersymmetric models, where the gluino-quark-squark vertex can be flavour violat- 
ing, the exchange of gluino and squarks in the loop gives contribution to the decay b ^ sj. 
Various combinations of the gluino-quark-squark vertex lead to |A(i?)| = |A(S')| = 1 mag- 
netic and chromomagnetic operators (of Oi-tjpe, Og-type) with an explicit factor a^, and 
to four-quark operators, with a factor a^. The complete effective Hamiltonian can then be 
split in two terms: 

^e// = + ^f// , (9) 

where 'H^jj- is the SM effective Hamiltonian in (^) and '^^fjj originates after integrating out 
squarks and gluinos. Note that 'mixed' diagrams, which contain, besides a W boson, also 
gluinos and squarks, give rise to as corrections to the Wilson coefficients in (at the 

matching scale). Such contributions can be omitted in a LO calculation, but they have to 
be taken into account at the NLO level. 

As far as the gluino-induced contribution to the decay amplitude 6 ^ 57 is concerned, 
the aim is to resum the following terms: 

(LO): a.(«.L)^, (iV = 0,l,...) 
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— ^ ^ — 

bn hi g 



FIG. 1. Diagram mediating the 6 — > ,57 decay through gluino exchange and contributing to the 
operator Ojb^g. A contribution to the primed operator Oy^g is obtained by exchanging L R. 

(NLO): «,«,(«,L)^, 

respectively at the leading and next-to-leading order. 

While 'Hljrjr is unambiguous, it is a matter of convention whether the factors, peculiar 
of the gluino exchange, should be put into the definition of operators or into the Wilson 
coefficients. In analogy to the decay b — > si^i~ discussed above, it is convenient to distribute 
the factors of as between operators and Wilson coefficients in such a way that the first 
two of the three steps in the program for the SM calculation also apply to the gluino- 
induced contribution. This implies one factor of al in the definition of the magnetic and 
chromomagnetic operators and a factor in the definition of the four-quark operators. 
With this convention, the matching calculation and the evolution down to the low scale 
/ife of the Wilson coefficients are organized exactly in the same way as in the SM. The 
anomalous-dimension matrix, indeed, has the canonical expansion in as and starts with a 
term proportional to a]. The last of the three steps in the program of the SM calculation 
requires now an obvious modification: the calculation of the matrix elements has to be 
performed at order as and a^ at the LO and NLO precision. With this organization of QCD 
corrections, the SM Hamiltonian Ti-^ff in eq. and the gluino-induced one Hfjj undergo 
separate renormalization, which facilitates all considerations. 

The effective Hamiltonian '^^f//, is further split into two parts: 

^!// = E QMOM + E E ci^{f^)oUf^) , (10) 

i i g 

where the index q runs over all light quarks q = u,d, c, s, b. The operators contributing to 
the first part are: 

• magnetic operators, with chirality violation coming from the 6-quark mass: 
07b,g = e g^sif^) m,(/i) (sa'^-'PRb) F,, , 0',,-^ = e gl{^x) m,(/i) {m^^'^PLb) F,, , 

0,,-g = gsifi) glilA m,{^) (sa^^T^PRb) G% , O^,,^ = (/x) mb(/i) (sa'^^T'^PLb) G% , 

(11) 

of dimension six, as the SM operators. A contribution to the magnetic operator 07b,g is 
shown in Fig. |l|. In this and the following diagrams, only the first in the series of possible 
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^ ® ^ 

bu g sl 

FIG. 2. Contribution to O^g^g from the insertion of the gluino mass and of a scalar mass term 
simultaneously violating chirality and flavour. A contribution to Ojg g is obtained through the 
interchange L <^ R. 




' ® ' ' ® ' 

bR g Sl bn g sl 

FIG. 3. Contributions to OYg^g from the insertion of the gluino mass and distinct chirality- and 
flavour-violating scalar mass terms. In the approximation nig = 0, the second diagram requires 
trilinear terms not linked to Yukawa couplings. The analogous contributions to Ojg g are obtained 
through the interchange L <^ R. 



insertions of chiral-flavour-violating scalar mass terms is drawn. This has the advantage of 
showing pictorially the correlation among supersymmetric sources of flavour violation and 
the generation of operators contributing to the effective Hamiltonian ([T0|). Nevertheless, the 
actual calculations presented in this paper are performed using squark mass eigenstates, i.e. 
resumming over all possible scalar mass insertions. 

• magnetic operators in which the chirality- violating parameter is the gluino mass rrig, 
included in the corresponding Wilson coefficients: 

Notice that these operators have dimension five, i.e. dimensionality lower than that of all 
remaining operators, of dimension six. Diagrams generating these operators are shown in 
Figs. and ||. 

• magnetic operators, with chirality violation signalled by the presence of the c-quark mass: 

Ojcrg = e gl{tj) rn.ifi) (sa'^'^PRb) F^, , O'j.^g = e gUii) m,{fi) (sa^^'PLb) F^, , 

0,c,g = gsifi) gUfi) (sa^'^T^PRb) G^^ , O'^^-g = g^ifi) g'M m,{fi) {sa^^T-PLh) G^, . 

(13) 
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FIG. 4. Penguin diagram contributing to the operators (|T 

The origin of these will become clear after discussing the second term in (|T0|) . This contains: 
• four-quark operators with vector Lorentz structure: 



Cl2,g = gt{^^){Sal^.PLhp) {qpl^Piqa) , C»12,g = 9M{salt^PRbp) [q^l^Phqa 

Ol,,g = gMirs^.PLb) (q^PRq) , OU,-g = gMis^.PRb) {q^Piq) , 

Ci4,g = gtif^){Salf,PLb/3) {qpYPRqa) , = g^^ (fi) {s a^f ^.Pub (3) {q^YPiqa 



(14) 



where colour indices are omitted for colour-singlet currents. They arise from box diagrams 
through the exchange of two gluinos and from penguin diagrams through the exchange of 
a gluino and a gluon. A typical penguin diagram is shown in Fig. ^. According to their 
Lorentz structure, these operators will be called hereafter vector four-quark operators. 

• four-quark operators with scalar and tensor Lorentz structure: 





= gtifx)isPRb) (qPnq), 


'-^15,5 


= gMisPLb) (qPLq), 




= 9t{^^){saPRbf3) (qfsPRqa) , 


•-^16,9 


= 9t{f^){saPLbi3) {q^Piqc) , 




= gtifx)isPRb) (qPLq), 




= gt{^i){sPLb) (qPRq), 


(nq 


= 9t{fJ'){saPRbf3) (qf^Piqa) , 


•-^18,9 


= 9t{f^){SaPLbi3) {qpPRqcc) , 


(nq 


= gt{fi){sa^uPRb) (qcr'^'PRq), 


'-^19,9 


= gt{^i){sa^M) (qcx'^-'PLq), 


(^q 


= g1{^j){sc(y^,uPRbp) {qpcr^" PRqa) , 


'^20,g 


= gl{^){sa(y^uPLbp) {qpcr^^PLqa 



(15) 

which are induced by box diagrams only and through the exchange of two gluinos. Examples 
of box diagrams are sketched in Figs. ^. In the following, the operators (|T5[) will be called 
scalar/tensor four-quark operators. Notice that, for different g's, Oii g-02Q^g are in general 
distinct sets of operators. 



The four-quark operators in ([T^) and (15) are formally of higher order in the strong 



coupling than the magnetic and chromomagnetic operators ([TT])- (0). As it will be explicitly 
shown in Sec. the scalar/tensor operators 0\^ g-02Qg mix at one loop into the magnetic 



and chromomagnetic operators. Given this fact, the necessity of including 0-jc,g and 0^c,g 
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in the operator basis becomes clear immediately: some of the operators Of^ g, Cgo^g with 
q = c mix into C7c,g and Csc.g- Such mixing terms can be calculated by considering the one- 
loop matrix elements (s7|Of_^|6) and {sg\Oig\b) {i = 15, ...,20), respectively. In principle, 
also operators like Oiu,g, 07d,g and O^s^g are induced in an analogous way. These operators, 
however, are weighted by m^j, irid and iris and are vanishing in the approximation used here: 
rriu = rrid = rris = 0. 

Due to these mixing effects, the scalar/tensor operators have to be included in a LO cal- 
culation for the decay amplitude. The remaining four-quark operators with vector structure 
Oil g^^iA g (and the corresponding primed operators) do not mix at one loop neither into 
the magnetic and chromomagnetic operators nor into the four-quark operators g-02Q g- 
Therefore, these vector four-quark operators become relevant only at the NLO precision. 

We end this section with a comment on the definition of the strong coupling constant 
used in the various steps of the calculation. In the full theory, which consists here of the SM 
and gluino-down-squark sectors of a supersymmetric model, all particles contribute to the 
running of this coupling, indicated by the symbol QsifJ')- In order to perform the matching 
with the effective theory, where only the five light quarks survive, all the heavy particles have 
to be decoupled. The strong coupling constant in this regime, indicated by (y's(/i), differs 
from Qsi/J') by logarithmic terms signalling the decoupling of the heavy particles: 

fi's(/i) = fi'.(/i) [1 + fi-^l/i) (decoupling log's)] . (16) 

At NLO precision, these decoupling terms have to be taken into account explicitly. At LO 
precision, however, ^s(/i) and QsifJ') can be identified and gsifJ') is here always understood to 
be the MS strong coupling at the renormalization scale /i, running with five flavours. 



III. SOURCES OF FLAVOUR VIOLATION 



Supersymmetric models contain all sources of flavour violation present in a Two Higgs 
Doublet Model of Type II, i.e. the vertices with a charged boson: un-dij-W^ and Un- 
dRj-H~^, URi-dij-H^ = 1,2,3). Once the electroweak symmetry is broken, a rotation 
in flavour space |^| 

D° = VdD, U° = VuU, D'° = UID\ W° = U*U% (17) 

of all matter superfields in the superpotential 

W = -Dr{hd),^QJHd + Ur{K),^Q°H^-fiHdH^, (18) 

brings fermions from the current eigenstate basis {ci^, m^, rf^, m^} to their mass eigenstate 
basis {dL,UL,dji,UR}: 

dl = VddL, ul = VuUL, d"ji = UddR, u°ji = UuUR, (19) 

and the scalar superpartners to the basis {D, U, D^, W^}. Through this rotation, the Yukawa 
matrices h^. and are reduced to their diagonal form and h^. 
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FIG. 5. Diagrams contributing to the operators (|T^). In the two upper diagrams, the quark 
(squark) q {q) can be of up- or down-type and the flavour violation on the lower squark line, 
not explicitly indicated, can be realized through a direct flavour-chiral transition (see Fig. ||) or 
through distinct chirality and flavour transitions (see Fig. In the lower diagrams, the down-type 
quark di is a h- or an s-quark if a single flavour violation is allowed in the squark lines. 

{hd)u = {UlhdVa)u = ^ , ih^u)u = {UlKV^)u = ^ . (20) 

Tree-level mixing terms among quarks of different generations are due to the misalignment 
of Vd and Vu, all the above vertices un-dLj-W^ and ULi-dRj-H~^ , ujn-dij-H^ = 
1,2,3) are weighted by the elements of the CKM matrix V = V^Vd- The supersymmetric 
counterpart of these vertices, un-Dj-W^ , ULi-D^*-H^ , URi-Dj-H~ , are also proportional 
to Vij in the limit of unbroken supersymmetry. 

To illustrate the sources of flavour violation that may be present in supersymmetric 
models in addition to those encoded in the CKM matrix, it is instructive to consider in 
detail the contributions to the squared-mass matrix of a squark of flavour /. The relation 
between off-diagonal terms in this squared-mass matrix and the type of operators inducing 
the decay b —>■ 57, will then become clear. Since present collider limits give indications that 
the squark masses are larger than those of the corresponding quarks, the largest entries in 
the squark mass matrices squared must come from the soft potential, directly linked to the 
mechanism of supersymmetry breaking. When restricted to the terms relevant to squark 
masses and quark-flavour transitions, the soft potential can be expressed in terms of the 
current eigenstates scalar flelds as: 

/i,,,i7,g^Df + A„,,i/„Q°[7;° + h.c.) . (21) 
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In ([21|), m^, m'^^, and m~ are hermitian matrices. The gluino g, a four-component Majo- 
rana spinor, is expressed in terms of the Weyl spinor A3 and has mass rrig = M3. Notice 
that, for the trihnear terms Ad,ijHdQ°Dj°, no proportionahty to the Yukawa couphngs is 
assumed. These trihnear scalar terms are left completely general and may also represent 
non-holomorphic ones, of the type A'^ijH*Q°Dj° discussed in |]36| . 

Thus, in the interaction basis (Qi, ^2°*^ Qs"*); often denoted also as (g£i, 

^^£25 ^£35 QRiy ?R25 Qrs)^ the squared-mass matrix for a squark of flavour / has the form 

2 _ ( /, LL + FfLL + Df LL (m^ ^ + Ff LR 

/ ~ / \ t 

V +FfRL m^j^jiji + FfRR + DfRR 

The term m j is , for both, up- and down-type squarks; m ^ is for a down-type 
squark and m ~ for an up-type squark. The off-diagonal 3x3 block matrix m j is A^Vd for 
a down squark, A*Vu for an up-type one. (The two vacuum expectation values are chosen to 
be real.) It should be stressed that, differently from '^/^ll and m'^f^RR, the off-diagonal 3x3 
matrix is not hermitian. In other words, it is A^^ij 7^ ^dji ^^^^ ^u,ij 7^ ^uji- 

The D-term contributions Dfn and Djrr to the squared-mass matrix (P^, 

/ = COS 213 M| (r^ - Q f sin^ ^h/) Us , (23) 

are diagonal in flavour space. 

The explicit form for the F-term contributions can be obtained from scalar quartic 
couplings arising from the superpotential (|18|): 

d vjDr{hihd)^^^ + viDr{hdhi)^^^"* - {i,v^Drh\^pY* + h.c.) + 

vlur {hiK) ^U"^ + vlur {Khi) ,^U'r"* - [^^VdUrhi,fi'^°* + h.c.) . (24) 
The rotation (|T7| ) reduces Ff ll and Ff rr to their diagonal form 

m^^D^A, mi,u:u,, nii.^rm, niifirut. 

as well as FfiR {Fj rl = fJ^j^) to 

-fi{md,i tan f])D*D^*, -^i{m^,^ cot f3)U*U^* . 

Therefore, once up- and down-quarks are brought to their mass eigenstate basis through 
the rotation (p!?]), the only sources of flavour violation in the squark sector arise from the 
off-diagonal terms in the soft mass matrices m ^ m ^ and m ^ [] Their origin, as 




^No new symbols are introduced to indicate the unknown matrices m?j rr?^ and m^- 
after the rotation ([l7|). Notice, however, that ^'^d equal before this rotation, are 



now related as = ^'^dhiX'^ ■ 
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their magnitude, is a model-dependent matter based on the interplay between the dynamics 
of flavour and that dictating the breaking of supersymmetry. In general, however, they give 
rise to large flavour-quark transitions at the loop level, through large couplings of gluinos 
to quarks and squarks belonging to different generations. 

One very drastic approach to this supersymmetric flavour problem is that of mSUGRA. 
In this model (or class of models) the soft potential (|2l|) is characterized at some high scale, 
typically a grand uniflcation scale, by the universality of the scalar masses: 

mpr.. = mr;.. = m ^ . . = m5ij\ (25) 

and the proportionality of the trilinear terms to the Yukawa couplings, through a universal 
parameter A: 

^d,ij = Ahd^ij ■^u,ij = Ahu,ij- (26) 

At this high scale, the only source of flavour violation is contained in the superpotential, 
indicating that the breaking of supersymmetry occurs at a scale where the dynamics of 
flavour has already taken place. 

An elegant solution to the flavour problem is obtained in GMSB models, in which the 
signal of supersymmetry breaking is transmitted to the visible sector of fields Q°, U°, D°, 
Hi,H2, etc., by flavour-blind gauge interactions. In these models, at the scale of supersym- 
metry breaking, all matrices in (p5D are diagonal, although different, and the common value 
of A in (|26|) is set to zero. 



In both mSUGRA and GMSB models, sources of flavour violation in the scalar sector are 
generated radiatively at the electroweak scale through the scalar quartic couplings propor- 
tional to Yukawa matrices. A simple inspection shows that intergenerational mixing terms 
due to only one type of Yukawa matrix, get eliminated by the rotation (|T^): no off-diagonal 
terms are therefore possible in m j in these models. On the contrary, flavour-violating 
terms are not rotated away in the m ^ ^ sector in which radiative contributions arise from 
quartic scalar couplings proportional to both matrices hd and hu. Being loop-induced, this 
source of flavour violation is, in general, small but it becomes non-negligible for large 
values of tan/3 [^]. By this reasoning it becomes clear that, while a contribution to the op- 
erator Orb,!} can arise from an off-diagonal term mixing the second- and third-generation left 
squarks {fn'^d, 11)23, as shown in Fig. |I], no contribution to Oi^f^g is possible in mSUGRA and 
GMSB models. The same holds for all other primed operators. These operators may never- 
theless acquire non-vanishing contributions in more general models, in which, for example, 
there exists an off-diagonal term (m^ rr)23- 

Also vanishing, in mSUGRA and GMSB models, is the contribution to the operator 
Ojg^g coming from a left-right mixing element (m^ lr)23- A contribution to this operator 
can however be induced, even in these models, by intergenerational mixing terms in j^^, 
(m^ 11)23, and the flavour-diagonal left-right term (m^ lr)33- ^^e mass-insertion formal- 
ism, often used for the calculation of supersymmetric contributions to FCNC processes |BS 



the first non- vanishing contribution to Ojg^g is then generated by the double insertion shown 
in the first diagram of Fig. ^ It will be shown later that, in generic supersymmetric models, 
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this contribution to C7§,g turns out to give the strongest constraint on {rn'^d^LLJ'iS: when 
reasonable values of (m^ lr)33 chosen. 

As advocated in the introduction, the aim of this paper is to provide a calculation as 
general as possible of the gluino contribution to the decay b — > 57, i.e. a calculation that 
applies to supersymmetric models with the most general soft terms. The QCD-corrected 
branching ratio for this decay can then be used to constrain the size of the off-diagonal ele- 
ments of the mass matrices tti^ j^j^, rri iti^ j^j^. Since different operators contribute 
to this decay, with different numerical impact on its rate, some of these flavour-violating 
terms may turn out to be poorly constrained. Thus, given the generality of such a calcu- 
lation, it is convenient to rely on the mass eigenstate formalism, which remains valid even 
when the inter generational mixing elements are large. The procedure used follows closely 
Refs. [^^. The diagonalization of the two 6x6 squark mass matrices squared Ad"^ and 
Aif^ yields the eigenvalues m^^ and {k = 1, ...,6). The corresponding mass eigenstates, 

Uk and dk {k = 1, 6) are related to the fields ulj, urj and dij, duj (j = 1, 3) as: 

Ul,R = ^UL,R U , dL,R = ^DL,R ^ , (27) 

where the four matrices r^2.,_R and F^i^ are 6x3 mixing matrices. The gluino-quark-squark 



vertices are explicitly given in Ref. [ p6[] . 



IV. WILSON COEFFICIENTS AT THE ELECTROWEAK SCALE 

At the matching scale fiw^ the non- vanishing Wilson coefficients for the SM operators in 
eqs. (IP and (|^) are, at leading order in a^: 

C2(/Uvy) = 1 

^^^^^'> ^ 24 (xt^""- 1)4 + ^^tw + ^^^i"' ^^^^t^ ~ ^^^^^'^ ^"^Xtn,) 

Csifiw) = g (j;^*"", 1)4 + ^^t^ ~ ^^^"^ ^ ^ ^ ^^^"^ Inxt^) , (28) 

with = m'f/M^. 

Among the coefficients arising from the virtual exchange of a gluino at the matching 
scale, the non-vanishing ones are |: 



C7,,sM = CiR) Y: (r^^^ T},'i) F,{x,dk) 

k=l dk 

CrsM = ^9^2 CiR) E ^ (r^^« T},'£) F,{x,dJ , (29) 

fc=l dfe 



^The linear combination Cjh^g{fj,w)C>7iy^g{fj,w) + C7g^g{fiw)C^7g,g{tJ'w) coincides with the expression 
5CT{iJ.w)Oj{fj,w) given in the literature (see e.g. ||2^ , |39|| ), where is the standard model operator. 



15 



in the case of magnetic operators and 

fc=i 4 

Css,^M=m~,-^j2^{^DR'^Di) [(C(i?)-iC(G))F4(x,,J-iC(G)F3(x,,J] , (30) 

in the case of chromomagnetic operators. The coefficients CYg^g{^w) and Csg,g{^w) of higher 
dimensionahty to compensate the lower dimensionahty of the corresponding operators. The 
ratios Xgd^. are now defined as Xgd^ = rri^/m'^j the Casimir factors C{R) and C{G) are 
respectively C{R) = 4/3 and C{G) = 3; and the functions -Fj(x), i = I, ...,4, are given in 
Appendix |^. The Wilson coefficients of the corresponding primed operators are obtained 
through the interchange F^^ij ^ ^^dl ^qs. ( p9|) and (|30D . The coefficients of the magnetic 
and chromomagnetic operators, proportional to the c-quark mass, vanish at the matching 
scale at lowest order in 

Compared to the SM, there is a larger number of magnetic and chromomagnetic operators 
with different chirality and dimensionality. The different chiralities are due to the fact that 
the gluino couples both to left- and right handed quarks and the associated squarks. In 
contrast, the W has only left-handed couplings and therefore right handed fields only arise 
if their masses are not neglected; usually only (chromo)magnetic operators with right-handed 
fe-quarks are included. Similarly, the occurence of (chromo)magnetic operators with differing 
dimensions can also be understood from the chirality structure of the gluino couplings. Some 
of the new operators differ from the SM (chromo) magnetic operators only by an additional 
factor (yf^. These were introduced as additional operators for practical reasons. 

Penguin diagrams mediated by the virtual exchange of a gluino and a gluon, yield non- 
vanishing coefficients only for the operators Oli g-Oi^ f 

ClgM = ^ E ^ (r^L n'l) mR) -\C{G)) FeixgdJ + iCiG)F,ixgdJ] 

k=l dk 

Ct2raM = -r^ E^(r^LrLL) [{G{R) -iG{G))F,{xgdJ + kC{G)F,{xgd,)] 
iOTT ^^-^ 

Cf-i^gif^w) = CfigifXw) 

Cl4,gi^^w) = Cl2,g(/^w) ; (31) 

as well as coefficients for the corresponding primed operators, 0\[g-0\'^g, which can be 
obtained from those in eq. (|3l|) by interchanging r)^^ <-> T^^i- These coefficients are actually 
independent of the quark label q. 

Box diagrams 0, with exchange of two virtual gluinos, yield the following contributions 



^Note that these diagrams are finite and all the manipulations needed to eliminate the charge 
conjugation matrices in the crossed topologies shown in Figs. @ are well defined. 
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to the coefficients Cfi g-Cf^^ g-. 



^i\^(M = T7^;i E {^{^'Ln'!){^tl^%\)[^G{x,,g,x,J + AF{x,,g,x,^^^^ 

-^"^ '""'9 k,h=l ^ 



+ ^.^4 (r^^r^^L) (r^r^^L) -20F(a;.,„x,,,)]| 

1 1 r 1 

ClS.gl/^w) = TTT^ "2 E 1 {^DL^*Dl) {^QR^Qr} [^G{Xd^g,Xq^g) — F{Xd^g, Xq^g)] 



+ ^idT:^ {^DL^Dl) {'^QR^Qr) G{Xd^g,Xqi^g) 



12 

1 1 ^ r 1 

Gl4,g{^^w) = E I e (-'^DLrDL) (J^QR^Qr) [G{Xd^g, Xq^^g) + 7F{Xd^^g, Xqi^g)] 

+ {^DL^h^l) {^I^R^Qr) G{Xd^:g,Xq^g)^ , (32) 

with the corresponding primed coefficients obtained through the interchange T^^)^ ^ ^^dr 
and Tqi^ ^ ^qr- Notice that the symbol 6qd is the Kronecker delta, equal to one when q is 
the down-quark and zero when g is a different quark. For q = d, also, the subscript Q in 
the two combinations (Tq^l^'ql) and (^q^r^'qr) has to be identified with D, typical of a 
down-type squark exchanged in the box diagram. The box-diagram functions G{x, y) and 
F{x, y) are explicitly listed in Appendix 0. 

The remaining coefficients Cl^g-C^Q^g-, in mass insertion language, are characterized by 
an odd number of L-R insertions in each squark line. In the mass-eigenstate basis used for 
squarks in this analysis, they are: 



91 ^1h9) 



1 I f 11 

Ci5,g(/^Ty) = Y^^^E^Iys ^dr^*dl) {^*ql^qr) Fi^d, 

g 

~ ^'J'^ 2 {^DR^*Dl} {^QL^Qr) Fi^dk9^^qh9) 

1 1 ^ ( 5 

Cte,gi^w) = Yq^^ ^ I 6 {^DR^*Dl) {^QL^Qr) F{^dk9-,Xqh9) 



"9 k,h=l 



+ ^dU^Dn^Dl) {^qI^q'r) H^d. 



k9i ^qh9) 
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1 1 ^ ( w 

Cu^gif^w) = H 1 {^DR^*Dl) {^h^R^QL) ^{Xd^g.Xq^g) 

'"'9 k,h=l ^ ^° 

~ '^S'^ 3 {^DR^*Dl) (J'qR^Ql) [^i^di:g,Xq^g) + 7F{Xd^:g^Xq^g)]^ 

1 1 ^ r 5 

Cis,g{f^w) = TZZ2Z^ X! 1 ~e ij''DR'^*DL) {JqR^Ql) G{Xd^g,Xg^g) 
'"'g k,h=l ^ " 

+ ^grf g {^DR^^J^l) {^QR^Ql) l^^i^dkgjXqhg) ~ ^i^dkgi^qhg)]^ 
1 1^(1 

Ci9-g{f^w) = 77-^— 2 J2 \ -o {^Dr'^*Dl) {^*q'l'^Qr) ^i^d^g^Xq^g) 
'"'9 k,h=l ^ ° 

+ ^id— (J'dR^d'l) {^QL^Qr) ^i^dkgjXqf^g)^ 

1 1 ^ r 3 

C20,g{f^w) = Vr-^~~2 \ o {^DR^*Dl) {^IJ'l^Qr) ^i^d^g^Xq^g) 

'"'g k,h=l ^ ° 

- 4 {^DR l) (rgl^ T%\) Fixd^g, xq,g) I . (33) 



The considerations made for the coefficients (|3^) hold also here: the corresponding primed 
coefficients are obtained through the interchanges T^dl ^ ^^dr ^ql ^ ^qr^ Sqd 
always vanishes, except for q = d. Under renormalization, the operators corresponding to 
the coefficients (0) mix with the magnetic and chromomagnetic operators in ( pT]) and (|13D 
by undergoing a chirality ffip proportional to m^. Therefore, only q = b and q = c can 
contribute to the decay 5 — »• 57 in the approximation of massless light quarks made here. 



V. WILSON COEFFICIENTS AT THE DECAY SCALE 



As already mentioned in Sec. fT|, the two terms T^Yff and 'Hlj-jr in the effective Hamiltonian 
(^) undergo separate renormalization. The anomalous-dimension matrix of the SM operators 
Oi-Os and the evolution of the corresponding Wilson coefficients to the decay scale fib are 
very well known and can be found in pF 



The evolution of the gluino- induced Wilson coefficients Ci^g from the matching scale fj,w 
down to the low-energy scale fib is described by the renormalization group equation: 

= Cj,g{fi) -fjirgifJ') ■ (34) 
The usual perturbative expansion for the initial conditions of the Wilson coefficients, 

Q,gM = cl^M + cl,M + , (35) 
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as well as for the elements of 'Jjigifi), 



= + + , (36) 



is possible thanks to the choice of including appropriate powers of gs{^^) into the definition 
of the operators Oi^g, as discussed in Sec. Since no NLO results are presented in this 
paper, the symbol 7ji,g(yu) will be used in the following to indicate the LO quantity 7jj^g(/x). 
Similarly the Wilson coefficients Ci^g will be indicating (7°^, as already understood in the 

361) run over all gluino-induced operators: 12 



previous sections. The indices i,j in (|35D and 
magnetic and chromomagnetic operators and 5 times (one for each flavour q) 20 four-quark 
operators. The anomalous- dimension matrix 7jj^^ is then a 112 x 112 matrix. It turns out, 
however, that primed and non-primed operators do not mix. This reduces the problem to 
the evaluation of two identical 56 x 56 matrices. 

Moreover, given their lower dimensionality, the dimension-five operators Ojg^g, Osg,g, and 
O'jgg, O'gg g, 6.0 uot mlx wlth d l mc u s lo U- s l X magnetic operators. The 4x4 submatrix for 
these operators is a block-diagonal matrix with 2x2 blocks. The block corresponding to 

07g,g, Osg,g is: 



18 
32 50 
~9 Y 



(37) 



and differs from the known mixing matrix of the SM operators Or and Og just by anomalous 
dimensions of the explicit mass mf, and of the coupling g'^ in the definition of the operators. 

In general, the structure of the remaining 54 x 54 matrix, corresponding to the four- 
quark operators O^g (z = 11, 20; q = u,d,c, s,b), magnetic operators On^g, Ojc^g, and 
the chromomagnetic operators Osb,g, Osc,g, is rather comphcated. The fact that in a LO 
calculation only the coefficients C7b,g and C7c,g (and corresponding primed coefficients) are 
needed at the low scale ^b, however, simplifies the analysis considerably. Among the four- 
quark operators, only those with scalar/tensor Lorentz structure, i.e. Ofg {i = 15, ...,20), 
mix into the magnetic and chromomagnetic operators at order as- The vector operators 
{O'lg {i = 11, 14)) on the other hand mix neither into the magnetic and chromomagnetic 
operators nor into the scalar/tensor four-quark operators. (The scalar/tensor operators, 
however, mix into the vector four-quark operators.) This implies that the presence of the 
four-quark operators with vector structure is completely irrelevant for the evolution of the 
coefficients of the magnetic operators. The observation that the scalar/tensor operators 
with the label q mix into Ojq^g and Osq^g, with the same q, together with the fact that 
scalar/tensor operators mix among themselves in a fiavour-diagonal way, further simplifies 
the situation. It is indeed possible to restrict the problem at the LO level to the calculation 
of two 8x8 matrices, i.e. the two matrices corresponding to the operators Ol^ g,OlQ g, 
Ci7,g, Ols,g, Ofg^g, Olo-g, Ojq-g, O^q-g, foT q = b aud q = c. 

For the case q = b, the result of such a calculation, in which the anomalous dimensions 
due to the explicit powers of the coupling as are again included, is: 



19 



44 

y 

-6 



16 





98 

y 




-48 







44 

y 

-6 









98 

y 





3 

1 

~2 




36 



-1 
7 
6 








-24 -56 



6 1^ 








1 

~3 
-1 





28 

y 

20 

y 

26 
32 




74 

y 



(38) 



The anomalous-dimension matrix corresponding to the case q = c differs from the previ- 
ous one in the submatrix responsible for mixing of the four-quark operators into the magnetic 
and chromomagnetic operators: 
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(39) 



Using the anomalous dimensions matrices (PT]), (^) and (PU]), the renormalization group 
equation (|34D can be solved by the standard procedure, described, for example, in Ref. pO| , 
using the Wilson coefficients Ci^g{iiw) given in Sec. fV|as initial conditions. The integration 
of (0) for Cjg^g and Csg^g yields the following expressions for these Wilson coefficients at 
the low scale fib- 



27 O / 25 27 \ 

25 



(40) 



Here and in the following, rj denotes the ratio as{fiw)/c(s{fi'b)- The low-scale Wilson co- 
efficients for the corresponding primed operators are obtained by replacing in (^0]) all the 
unprimed coefficients with primed ones. The same holds for the following coefficients. 
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The Wilson coefficients of the dimension-six operators Cii),g and C^h,g are at low scale: 



C7b,g{lJ'b) = ^723 Crb,g{l^w) + - (r^M - 7723) Csb,g{l^w) + RibrgifJ'b) , 

37 

Csb,gilJ'b) = C8b,g{lJ'w) + -Rsb.gl/^fc) ■ 



(41) 



The remainder functions R7b,g{fJ'b) and Rsb,g{,IJ'b) are given in Appendix ^. They turn out to 
be numerically very small with respect to the other terms on the right-hand sides of (|4l|) . No- 
tice that, in the approximation R7b,g{fJ'b) = R?,b,g{l^b) = 0, the low-scale coefficients Cn^gifJ'b) 
and Csb^gifJ'b) are simply obtained through the integration of (p^ ) with the anomalous di- 
mension matrix 'jji^g reduced to the 2x2 block of ( PS|) corresponding to the operators 07fe,g 
and 08b,g- 

Finally, the coefficients Cjc,g{,l^b) and C^cgil^b) formally have the same expression as 
C7b,g{fJ'b) and Csb,g{fJ'b), when the indices 7b and 86 are replaced by 7c and 8c. Also in this 
case, the functions RrcgifJ^b) and R^^gil^b)^ listed in Appendix ^ are numerically small. In 
the approximation RrcgifJ^b) = R8c,g{fJ'b) = 0, the coefficients CjcgifJ'b) and CgcgifJ^b) vanish 
identically, since the corresponding Wilson coefficients at the matching scale are vanishing. 



VI. BRANCHING RATIO 



The branching ratio BR(i? —>■ Xs'y) can be expressed as 

J- SL 



(42) 



where BR^^ = (10.49±0.46)% is the measured semileptonic branching ratio. To the relevant 
order in a^, the semileptonic decay width is given by: 



™5 IT/ 12 / 

rrif^Lrpl Vcb\ 
1927r3 




(43) 



where the phase-space function g{z) is g{z) 
width for b S'y reads: 



+ 8z^ - z^ 



T{b 57) 



mlGl\VttV;fa 
327r4 



+ 



C7 



Vlz^ log z. The decay 



(44) 



where C*7 and C'-j can be expressed in terms of the SM and gluino-induced Wilson coefficients 
evolved down to the decay scale as: 

16y27r3a3(/i6) 



16v^7r3a,(/i5) 



GpVtbV^ 



ts 



1 

Cjb.gifJ'b) H CYgg{fib) H C-jcA^-b) 

rrih 



(45) 



Notice that, at the leading logarithmic level, it is not possible to distinguish between the 
pole masses rrib and rric from the corresponding running quantities at the scale m?, or rric. In 
the following, these mass parameters are always treated as pole masses. 
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VII. NUMERICAL RESULTS 



Numerical predictions for the QCD-corrected branching ratio BR(5 Xsj) induced by 



gluino-squark exchange can be obtained from eqs. (^)-(|45|). To show these results, it is 
convenient to select one possible source of flavour violation in the squark sector at a time 
and assume that all the remaining ones are vanishing. 

Following Ref. [^, all diagonal entries in ^j^, and m\ are set to be equal 



and their common value is denoted by m?. The branching ratio can then be studied as a 
function of only one off-diagonal element in and normalized to m?, i.e. as a 

function of one of the elements 

7^ j) (46) 



JLL,ij — 2 ' "RR,ij — 2" 



and/or of one diagonal or off-diagonal element of the 3x3 matrices 

LRi ^d RL SS^-IH 



normalized to m?: 



(m 



\LR)ij 



'LR,ij - —2 , dRL,ij - —2 • (4g 



The corresponding off-diagonal entries in the up-squark mass matrix squared, relevant for 
the contributions coming from the gluino-induced four-quark operators (14) and ([T5|) are 



set to be equal to those in the down-squark mass matrix squared. Among the four-quark 
operators, only the scalar/tensor operators (|15]) contribute to BR{B — >• Xg'j), at the LO 
order in QCD. Their effect is negligible and the above restriction is not likely to produce an 
unnatural reduction of their contribution. Indeed, due to their proportionality to T']^^T*j^[, 
the operators Of g {i = 15, 20) are generated always together with O^g^g and Ogg^g. As will 
be discussed later, the latter are the numerically important operators and the corrections 
induced e.g. by Csg.g on the Wilson coefficent C^g^g^fib) completely overshadow the effect of 
the four-quark operators Olg. These induce corrections of the Wilson coefficient Cji-{fib) 
of the numerically less relevant operator Oji- that are generically suppressed by a factor 
{mb/m~) at the amplitude level. Analogously, the primed scalar /tensor operators O^'g {i = 
15, ...,20) are also expected to have a very small impact on the decay amplitude. The 
vector four-quark operators, on the other hand, can be generated without the simultaneous 
generation of O^g^g and Csg.g and no suppression factor {nib/m-) is present in this case. 
Therefore, the vector four-quark operators, although entering at NLO only, are in general 
expected to have a larger impact on the decay amplitude than the scalar/tensor four-quark 
operators. In the context of a NLO analysis, one should actually check if the assumption 
of equal off-diagonal entries in the up- and down-squark mass matrices squared is not an 
oversimplification, affecting the generality of the numerical results. 

As for the remaining entries in the squark mass matrices squared, the Z)-terms are 
calculated using Mz = 91.18 GeV, siia^ 9w = 0.2316, and tan/3 = 2; the F-terms -F/ ll and 
Ff RRi using nib = 3 GeV and = 175 GeV, in the approximation of vanishing lighter quark 
masses, whereas Fj lr = Ff = is assumed. It is obvious that all the information gained 
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FIG. 6. Gluino-induced branching ratio BR(5 — > ^^7) as a function of x = m?/m?, ob- 
tained when the only source of flavour violation is 6lr,23 (see text), fixed to the value 0.01, for 
rriq = 500 GeV. The solid line shows the branching ratio at the LO in QCD, for fi^, = 4.8 GeV 
and nw = ^w] the two dotted lines indicate the range of variation of the branching ratio when 
fj-h spans the interval 2.4-9.6 GeV. Also shown are the values of BR(5 Xgj) when no QCD 
corrections are included and the explicit factor as{^) in the gluino-induced operators is evaluated 
at 4.8 GeV (dashed line) or at Mw (dot-dashed line). 



V3 

o 



10 



X 

A 
I 

PQ 
PQ 



1 - 



"T — I — r 



"1 — I — I — I — I — I — I — I — I — I — I — I — I — r 

gluino only 



J I L 



J I I L 



J I I L 



J I I L 



0.5 



1.5 



FIG. 7. Same as in Fig. when only 6ll,23 is non-vanishing and fixed to the value 0.5. 
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through the numerical evaluation of BR(S Xs'~f) on the size of (m^ lh)33 "^^^ t)e extended 
to the combination (m^ + FfLji)^^ and ((^^.2,/?)^ + -F/i?L)33 in realistic cases, in which 
/X 7^ 0. For the diagonal entry m?, the value rrig = 500 GeV is in general used. Moreover, 
it is imposed that the eigenvalues of the two 6x6 up- and down-squark mass matrices are 
larger than 150 GeV for all values of the 5-ratios scanned. The value of 150 GeV is here 
taken as an average model-independent lower limit on squark masses, which can be inferred 
from direct searches of squarks at hadron colliders. 

Finally, the remaining parameter needed to determine the branching ratio is: 

2 

X = ^ , 48 

where rrig is the gluino mass. 

In the following, the SM contribution to BR{B — > Xg'j) is, in general, added to the gluino 
contribution: possible constraints on the flavour- violating sources in the squark sector should 
be extracted, keeping into account that the SM contribution already successfully saturates 
the experimental result for this branching ratio [|l^,|l^. As already stressed in Sec. |, this 



analysis applies to particular directions of the supersymmetric parameter space, in which 
charged Higgs, chargino and neutralino contributions can be safely neglected with respect 
to the gluino and SM contributions. Moreover, it should also be mentioned that the bounds 
discussed in this section on Sll^23, Srr,23, Slr,23, and Srl^23, obtained in these particular 
directions of parameter space, have to be understood in an indicative sense, since they are 
extracted ignoring the error of the theoretical calculation. 

It is useful to isolate the gluino contribution when illustrating the impact of the LO 
QCD corrections on the gluino-induced Hamiltonian. In Figs. ^ and 0, indicated by solid 
lines, are shown the values of the QCD-corrected branching ratio obtained, respectively, 
when only Slr,23 and 5ll,23 are non-vanishing. Their values are fixed in the two Figures as 
follows: (5iR,23 = 0.01 and 5ll,23 = 0.5. The branching ratio is plotted as a function of x, 
i.e. as a function of the gluino mass, for a given value of rriq, rriq = 500 GeV. Also shown is 
the range of variation of the branching ratio, delimited by dotted lines, obtained when the 
low-energy scale fib spans the interval 2.4-9.6 GeV. The matching scale fiw is here fixed to 
Mw- As can be seen, the theoretical estimate of BR(5 — > Xsj) is still largely uncertain 
(~ ±25%). An extraction of bounds on Sll,23 and Slr,23 more precise than just an order 
of magnitude would require, therefore, the inclusion of NLO QCD corrections. It should be 
noticed, however, that the inclusion of corrections at the LO has already removed the large 
ambiguity on the value to be assigned to the factor as{fi) in the gluino-induced magnetic 
operators ([TT|)-(p!3D. Before adding QCD corrections, the scale in this factor can assume all 
values from fib to nw The corresponding values for BR(5 — ^ Xs'j) for the two extreme 
choices of /x are indicated in Figs. ^ and |^ by the dot-dashed lines (/i = Myy) and the dashed 
lines {fi = 4.8 GeV): the branching ratio is virtually unknown! The choice /i = M^y gives 
values for the non-QCD-corrected BR.{B Xg'y) relatively close to the band obtained for 
the QCD-corrected result, in the case shown in Fig. 0, when only 6ll,23 is non-vanishing. 
Finding a corresponding value of fi that minimizes the QCD corrections in the case studied 
in Fig. 0, when only 61^23 is different from zero, depends strongly on the value of x. 
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FIG. 8. Dependence of the QCD-corrected branching ratio BR(i? ^sj), obtained from the 
SM and gluino contributions, on the parameter (5l_r,23) when (m^ lr)'23 only non- vanishing 

off-diagonal element in the down-squark mass matrix squared. The branching ratio is shown for 
different values of x = m|/m|, with nig = 500 GeV: 0.3 (short-dashed line), 0.5 (long-dashed line), 
1 (solid line), and 2 (dot-dashed line). Low and matching scales are: /ife = 4.8 GeV and iiw = M^y. 
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FIG. 9. Same as in Fig. |8|, when 6jiL,23 is the only source of flavour violation for the gluino 
contribution. The parameter x is fixed to: 0.3 (short-dashed line), 0.5 (long-dashed line), 1 (solid 
line), 2 (dot-dashed line). 
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The results in Figs. ^ and |^ also show that the operator On^g gives much smaller con- 
tributions to BR(i? — > Xg'-f) than the operator O^g^g. Indeed, the branching ratio obtained 
through OYh,g only is typically suppressed by a factor {mh/rrigY, with respect to that obtained 
from O'Tg^g, if similar values of 5ll,23 and Slji^23 are chosen. Analogous considerations hold 
for O'-^jjg and O'-^g g. The elements Slr,23 and Srl,23 are therefore expected to be the flavour- 
violating parameters most efficiently constrained by the measurement of BR(i? Xg'y). 

In Fig. ^, the dependence of BR(i? — > Xs'j) is shown as a function of 6lr,23 when this 
is the only flavour-violating source. The two horizontal lines correspond to the minimum 
and maximum values, 2 x 10"'^ and 4.5 x 10~^, allowed by the CLEO measurement. The 
branching ratio is obtained by adding the SM and the gluino contribution calculated for 
different choices of x, and a fixed value of mf rrig = 500 GeV. The values of the gluino 
mass corresponding to the choices x = 0.3, 0.5, 1, 2 are: rrig = 274, 354, 500, 707 GeV. The 
branching ratio is plotted in this Figure for fixed values of the two scales: fib = 4.8 GeV 
and fiw = Mw- The gluino contribution interferes constructively with the SM for negative 
values of 5l_r,23; which are then more sharply constrained than the positive values. Overall, 
this parameter cannot exceed the per cent level. No interference with the SM is present 
when 6rl,23 is the only source of flavour violation, as shown in Fig. ^ The results obtained 
for BR(i? Xs'~f) are then symmetric around 5rl,23 = and the constraints on |5/jl.23| are 
upper bounds on its absolute value: there are no small values of 5_rx,23 for which the total 
branching ratio falls off the band allowed by the CLEO measurement. 

Much weaker is the dependence of BR{B —> Xg'j) on 6ll,23 if ("^d ll)23 is the only off- 
diagonal element in the down squark mass matrix squared. This dependence is illustrated 



in Fig. |10| for different choices of x and rUg = 500 GeV. The gluino-squark loop generates 
in this case only the dimension-six operator 0^,13 and the gluino contribution interferes 
constructively with the SM contribution for positive Sll,23- Notice that the mass insertion 
approximation, given the large values of Sll,23 allowed by the experimental measurement, 
cannot be used in this case to obtain a reliable estimate of BR(S — > Xg'y), whereas it is 
an excellent approximation of the complete calculation in the cases shown in Figs. |] and |^. 
For completeness, also the case in which the only off-diagonal element in the down-squark 
mass matrix squared is in the right-right sector, (m^ rr)23 7^ 0, is shown in Fig. |Tl|. The 
inclusive branching ratio, plotted versus the relevant parameter 6rr^23, is now obtained from 
the incoherent sum of the SM and gluino contributions and shows conspicuous deviation 
from the SM result only for very large values of 5r_r,23- 

As already observed, among the operators C7b,g and Ojg^g, the second one has the stronger 
impact on BR(i? — > Xg'j). It is then legitimate to question whether O-^g^g may not provide 
a stronger constraint on 6ll,23- Since Ojg^g requires a chirahty flip within the loop, then at 
least an additional off-diagonal element different from zero is needed in the left-left sector of 
the down-squark mass matrix squared. Indeed, the flavour-conserving left-right mixing term 
[m^ l_r)33 together with (m^ ll)23 can also generate the operator C?^,^; see the first diagram 



in Fig. 0. The corresponding branching ratio is shown in Fig. |T2|, as a function of 611^23 for 
different choices of Slr,33- The value of the diagonal entries in the squark mass matrix is 
rriq = 500 GeV and m~ is determined by the choice x = 0.3. As in the previous plots, low 
and matching scales are fixed as /i;, = 4.8 GeV and fiw = Mw- Both parameters Slr,33 and 
Sll,23 are chosen to be positive. The solid line in this Figure, obtained for 5^^,33 = 0, then 



26 



h 6 



A 
I 

OQ 

PQ 







1 1 1 1 1 1 1 1 1 1 1 

- SM + gluino 


1 1 1 1 1 1 1 1/ ^ 

I r 

11 
/ /- 














^T'" : 



-1 



-0.5 







0.5 



LL.23 



FIG. 10. Same as in Fig. |8|, when 5ll^23 is the only source of flavour violation for the gluino 
contribution. The different lines correspond to: x = 0.3 (short-dashed line), 0.5 (long-dashed line), 
1 (solid line), 2 (dot-dashed line). 
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FIG. 11. Same as in Fig. ^, when b^.R^iz is the only source of flavour violation for the gluino 
contribution. The values of x corresponding to the different lines are: 0.3 (short-dashed line), 0.5 
(long-dashed line), 1 (solid line), 2 (dot-dashed line). 
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FIG. 12. BR(S Xg'y) vs. Sll,23, when Sll,23 and Slr,33 are the only sources of chiral-flavour 
violation. The dependence on 6ll,23 is shown for different values of Slr,33- (solid line), 0.006 
(short-dashed line), 0.01 (dot-dashed line), 0.1 (long-dashed line). The value of x = m^/rn^ is fixed 
to 0.3 and niq to 500 GeV. 



coincides with the short-dashed hne in Fig. 10. The SM value of the branching ratio, at the 



LO in QCD, is the value at which all curves meet for 5ll^23 = 0. The short-dashed line is 
obtained for (m^ Li?)33 — ''^g^b; which corresponds to a relatively large trilinear coupling 
in models in which the trilinear term in the soft potential is proportional to the Yukawa 
couplings. The corresponding maximally allowed value of 5ll,23 already is, in this case, 
considerably smaller than that obtained when only the operator 0-jh,g is present. Larger 
values of 5lr,33 obviously induce even more stringent constraints on 5ll,23- 

Two obvious lessons can be learned out of this analysis. First, in directions of the 
supersymmetric parameter space in which other contributions to BR(5 — »• ^^7) cannot be 
neglected, some of the constraints derived here may be invalidated by possible interferences 
among different contributions. An illustration of this is provided by the comparison of 
the bounds imposed by BR(i? — > Xg-^) on 5lr^23 and 5rl^23i which are different precisely 
because contributions from SM-gluino interferences are possible in one case, but not in the 
other. The second lesson stems from the observation that different operators contributing to 
BR(i? Xs"^) have very different numerical relevance. Because of this, it is not necessarily 
true that the strongest constraint on a chiral-flavour-violating sfermion mass term can be 
derived from the operator that is generated by it in the most straightforward way. Therefore, 
one cannot but end this section by stressing again the importance of analyses as complete 
as possible, when attempting to use the b —>■ sj decay as a model-building tool, constraining 
the soft supersymmetry-breaking terms. 
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VIII. SUMMARY 



Gluino-mediated contributions to FCNC processes are useful probes of chiral-flavour- 
violating soft breaking terms. They are in general cleaner than chargino contributions, 
which are sensitive also to the CKM matrix, responsible for flavour violation in the SM and 
in 2HDMs. Since they come with a coupling cts, they are usually rather large. Whether 
they are indeed much larger than chargino contributions is a model-dependent issue. 

The presence of the coupling ag makes these contributions also particularly interesting for 
FCNC processes in which QCD corrections play as important a role as the purely electroweak 
contributions. Exemplary among these processes is the decay h — * S7. A specific analysis of 
the implementation of QCD corrections for the gluino contribution to this decay is required. 
This paper is devoted to precisely this issue: it shows how to QCD-correct the gluino 
contribution to the decay & — > 57, using the formalism of effective Hamiltonian. 

It is shown here that, contrary to the common belief, gluino contributions require 
an enlargement of the standard basis of operators needed to describe h-s transitions 
in the SM and 2HDMs. In the SM, the calculation at the LO in QCD includes all 
terms of type {as\og{Mw /f^b))^ ■, whereas the calculation at the NLO resums all terms 
as^Oislog^My^/ml))^ . The program of implementation of QCD corrections in the SM re- 
quires that at each order in QCD, e.g. LO or NLO, the anomalous-dimension matrix of the 
SM operators is calculated at a higher order in ag than the matching conditions and the 
matrix elements. This is because in the SM, at a certain order in QCD, no other opera- 
tor can mix into the magnetic operator (e/(167r^)) {sa^^PRh) F^^, without the exchange of a 
virtual gluon. The situation is different in the case of gluino contributions. Gluino-induced 
magnetic operators acquire corrections as in the SM, when an additional virtual gluon is 
exchanged. Moreover, as in the SM, also gluino-induced chromomagnetic operators mix into 
the gluino-induced magnetic ones after the on-shell gluon is connected to a quark line and 
an additional photon is radiated. Both operators get first non-vanishing contributions at the 
matching scale at order a^, and give QCD-corrected contributions of type «^ log(M^/m^). 
Gluino-induced four-quark operators with first non- vanishing contributions at the matching 
scale of C(a^), however, can mix into the magnetic operators through the connection of 
two of the external quark lines and the emission of an on-shell photon, giving therefore also 
corrections of type log(M^/m5). As not all logarithms are due to gluon exchange, their 
systematic resummation is more involved as in the SM. 

A solution to this problem has been proposed in this paper. The couplings and 
al intrinsically connected with the gluino exchange are respectively factorized out in the 
definition of magnetic and chromomagnetic operators and of operators originating from 
box diagrams. With this definition, all gluino-induced operators are distinguished from 
the standard set of operators in the effective Hamiltonian induced by SM and 2HDMs. In 
particular, the magnetic operator O-j = (e/(167r^)) (so-^'^PRfe) F^^, is now distinct from 
the gluino-induced one Ojh^g = eglwib {sa^^Puh) F^y. This in turn has to be distinguished 
from the lower dimensionality operator Ofg^g = egl {sa^'^PRb) F^i,, induced at the matching 
scale by a loop diagram in which chiral-fiavour violation is provided, for example, by the 
insertion of a left-right mass term in the squark propagator and the insertion of a gluino 
mass in the gluino propagator. Completely new are the four-quark operators, such as ^15,3 = 
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gj{fj,){sPLb) (qPLq), with an explicit factor gff. In total, the inclusion of gluino contributions 
requires 56 new operators and other additional 56 with opposite chirality. 

With the above definition of gluino-induced operators, an important goal is achieved. 
The first non- vanishing contribution to the gluino-induced Wilson coefficients is of 0{a^^). 
Moreover, the anomalous dimension matrix starts at 0{as)- Consequently, the integration 
of the renormalization group equation yields, at first non- vanishing order, terms all of type 
{aslog{Mw /rrib))^ . The analogy with the LO SM contributions is now clear. It is this first 
non-vanishing order that is classified as LO gluino contributions. Thus, gluino exchanges 
induce terms of type as{as\og{Mw /rrii,))'^ at the LO in QCD, to be compared to the LO 
SM contributions of type Gf(q;s log(M^/m;,))^. The generalization to the NLO is obvious: 
it will yield contributions al{as \og{Mw / rrib))^ for the gluino-induced operators, versus the 
contributions G FO£s{o£slog{Mw /rrib))^ coming from the SM set of operators. 

A complete LO analysis for the branching ratio of the inclusive decay B — > Xg^ coming 
from SM and gluino-induced contributions is presented in this paper. The full anomalous- 
dimension matrix for gluino-induced operators is calculated and a simple expression for the 
branching ratio is given. The gluino-induced Wilson coefficients are also listed. They are 
obtained from the evaluation of one-loop diagrams mediated by the exchange of gluino and 
squarks. The mass eigenstate formalism is adopted as the most suitable for supersymmetric 
models with different sources of flavour violation and with a priori large flavour-violating 
mass terms. 

A numerical analysis for the inclusive branching ratio BR(S Xg'j) due to SM and 
gluino-induced contributions is presented. The QCD corrections to the gluino-induced con- 
tributions are found to be even more crucial than in the SM case. The non-corrected contri- 
butions to the inclusive decay B — > Xsj, in fact, suffer from a severe source of uncertainty 
that has no counterpart in the SM. At the zeroth order in QCD, there is no prescription to 
fix the scale of the overall factor cc^ in the final expression of the branching ratio, intrinsically 
due to gluino exchanges: it can range from the matching scale ~ Mw to the low-scale ~ m?,. 
Once QCD corrections are added, the bulk of this ambiguity is removed: this factor of 
has to be evaluated at a low scale of 0{mh), although the exact value of this scale remains 
unknown. A similar uncertainty is due to the fact that the matching scale is only known 
to be of 0{Mw)- Thus, the LO branching ratio still suffers from matching- and low-scale 
uncertainties similar in size to those in the SM results. 

Finally, we conclude by recalling that this analysis is valid in particular directions of the 
supersymmetric parameter space, in which charged Higgs, chargino and neutralino contribu- 
tions can be neglected. In spite of the still large theoretical error, it provides bounds on the 
different sources of fiavour violation that are present in these directions of parameter space. 
Further studies are called for to include NLO contributions as well as all the remaining 
supersymmetric contributions. 
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APPENDIX A: FUNCTIONS 



Listed below are the loop functions appearing in the coefficients ( p9D and (|30|): 
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those originated by the calculation of penguin diagrams (see coefficients (pT])): 
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APPENDIX B: WILSON COEFFICIENT REMAINDERS 



The effect of the four-quark 
relative to the magnetic and chroi 
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where the factors di-d^ are given by 
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di = rj 23 ; d^ = r] 23 ; rfQ = /y23. (B5) 



Notice that in (|B1D-([B^ there is no dependence on CiY g{fiw) and C'j'g ^(yUvi/), as there is no 
dependence on Cijg{fiw) and Cigg{iiw) in (P^)- ([B4D . By inspecting the two anomalous- 
dimension matrices in eqs. ( PBD and it is easy to see that the two operators Ofjg, Of^ g, 
do not mix with the remaining ones Ol^ g, Of^ g, Ofg g, C^cs' ^7g,5' ^Sq,g in either of the 
two cases, q = b and q = c. 
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